Laser generated wake fields as a new diagnostic tool for magnetized plasmas 
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In the presence of an external magnetic field, wake fields generated by a short laser pulse can 
propagate out of the plasma, and thereby provide information about the electron density profile. 
A method for reconstructing the density profile from a measured wake field spectrum is proposed 
and a numerical example is given. Finally, we compare our proposal with existing plasma density 
diagnostic tecniques. 



As is well-known, a short laser pulse propagating in 
an underdense unmagnetized plasma can excite a wake 
field of plasma oscillations This has interesting ap- 
plications to plasma based particle accelerators and is 
naturally of importance for the general understanding of 
laser-plasma interactions. If the plasma is magnetized 
and the external magnetic field non-parallel to the di- 
rection of propagation of the laser pulse, the wake field 
becomes partially electromagnetic and thereby obtains a 
nonzero group velocity 

In the present paper we study wake field generation 
in an inhomogeneous magnetized plasma. In particular 
we address the question to what extent the wake field 
can propagate out of the plasma, and thereby provide 
information about the plasma parameters. The most in- 
teresting case is that of a strongly magnetized plasma 
(i.e. when the electron cyclotron frequency is larger than 
the plasma frequency), for which almost all of the wake 
field energy - except that generated in a narrow low den- 
sity region - may propagate out of the plasma. Since 
the wake field also in the magnetized case initially has 
the frequency equal to the local plasma frequency, this 
provides a way of extracting information about the elec- 
tron density profile. It turns out that even though wave 
overtaking - for example when a higher frequency part 
of the wake field passes a lower frequency part - may 
occur, the density profile can still be reconstructed by 
integrating the ray equations of geometric optics back- 
wards. This is more straightforward for a static back- 
ground density profile, and a numerical example is pro- 
vided for this case, where the predicted spectrum of the 
wake field corresponding to an assumed density profile is 
shown, and a reconstructed profile is calculated. Finally, 
the method is compared with two of the existing plasma 
density diagnostic techniques 0, namely interferometry 
and refiectometry. 

We consider a high frequency laser pulse with fre- 
quency u>h propagating in a cold, weakly inhomogeneous 
magnetized plasma. We assume the ordering loh 3> 
Lo p ,uj c , where ui p and ui c = \q~Bo\/m are the plasma and 
electron cyclotron frequency respectively, q and to are 
the electron charge and mass, and Bo = Box is the (con- 
stant) external magnetic field. We let the laser pulse 
propagate perpendicularly to the external magnetic field. 



The ponderomotive force of the laser pulse will gener- 
ate a "low frequency" wake field mode (which is the low 
frequency branch of the extraordinary mode, or plasma 
oscillations modified by the magnetic field, depending on 
the choice of terminology) during its path through the 
plasma. The generation mechanism is most efficient if the 
pulse has a duration of the order of the inverse plasma 
frequency or shorter. In principle, the laser pulse will 
broaden due to ordinary dispersion, decrease its energy 
and frequency due to the interaction with the wake field, 
etc. These and other effects have been considered in ho- 
mogeneous plasmas by for example Rcf . . We will focus 
on the spectral properties of the wake field, however, and 
for this purpose it turns out that we can forget about the 
details of the laser pulse. Basically the effect of the laser 
field is to provide a well localized ponderomotive source 
term in the governing equations for the wake field, trav- 
elling with almost the speed of light in vacuum. 

The wake field quantities are denoted by index L. We 
introduce the corresponding vector and scalar potentials 
At,(z, t) and 4>l(z, t), using Coulomb gauge, and the elec- 
tron density is written n — no(z, i) + nr,(2, t), where no is 
the unperturbed density, assumed to vary on longer space 
and time scales than the wake field. Furthermore, the 
electron velocity v is divided into its high- and low fre- 
quency part, and we denote the low frequency contribu- 
tions perpendicular and parallel to the direction of prop- 
agation with Vii and vt. z respectively. The pondero- 
motive force of the laser pulse induces longitudinal wake 
field motion, which couple to motion in the y -direction 
through the Lorentz-force, but there is no wake field mo- 
tion in the direction of the external magnetic field, and 
accordingly we put Al = Ar,y and V£j_ = t>Lj_y. Lin- 
earizing in the low frequency variables, and neglecting 
derivatives acting on no, we obtain the following set of 
equations governing the wake field generation 
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Next we write this on matrix form Au = b, where u = 
(4>L,nL,VL Z ,VL±, Al), A is an operator matrix and b is 
the (ponderomotive) source vector. Since derivatives on 
no are neglected, we can simplify the system (|l])-(||) by a 
simple matrix inversion u =A _1 b treating the operators 
in A like constant coefficients. We introduce the group 
velocity of the laser pulse v g n = dkto and eliminate dif- 
ferential operators in the denominators everywhere by 
"multiplication". By using <9 t « ~v g nd z » — cd z for 
derivatives acting on the source vector, and focusing on 
the first component of u, we finally obtain after one space 
and one time integration 
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where u> 2 — lu 2 +u> 2 and we recognize the operator acting 
on <pi as the wave operator for the extraordinary mode. 
Although derivatives on no have been neglected in the 
derivation of Eq. (||) , it should be emphasized that the 
equations for wake field generation is more complicated 
here than for a homogeneous plasma. This is because in 
the inhomogeneous case, the evolution of the wake field 
is in general not quasi-static in a frame moving with the 
group velocity of the laser pulse, and thus dt ~ —v g iid z w 
— cd z does not hold for derivatives acting on the wake 
field. 

It is useful to divide the study of the wake field prop- 
erties into its excitation and its propagation phase. The 
excitation of one additional wavelength of the wake field 
takes place during a distance of the order of 2 , kc/uj p , and 
- as a basic assumption of ours - the variations of no is 
negligible on this length scale. Thus as far as the excita- 
tion process is concerned, the plasma can essentially be 
treated as homogeneous. The solution for the wake field 
can thus be obtained from previous authors 0. Chang- 
ing to comoving coordinates £ = z — v g nt, t = t the 
result for the wake field potential can be written 



<Pl = </>Losin[/c p (£ - f )] 
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p — LJ p /V g H, CO 



where k p is the wake field wavenumber k 
is the (constant) position of the (short) laser pulse, and 
4>lo = (qUp/mVg) J_ 00 | Ah 2 <i£. The important result 
here, for our purposes, is the determination of the initial 
value of the wake field wave number k p = uj p /v g H, which 
corresponds to an initial frequency u) p (in the laboratory 
frame), that will vary with the position of generation. 

Next we make the ansatz of geometric optics. Inspect- 
ing the wave operator in Eq. (^|), we note that there is 
a resonance at J 1 = ui 2 = u 2 + uj 2 and cut-offs at ui = 
\[-w e + K 2 + 4o$V2] and UR = i [uJc + (w a + 4^2)1/2], 
The dispersion relation has two positive roots. One 
branch is valid for uj > lor and one for u>l < u> < uJh- 
As the wake field is generated with the local plasma 



frequency u) p , the wake field must belong to the latter 
branch. Therefore we write the dispersion relation from 
now on as 



= W(k, Z, t)={ X - ^X 2 -^c2fc2_ CJ 4| 
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where x = ' 



-c 2 k 2 )/2. Due to the cut-off loi and 



the resonance u>h, only parts of the wake field satisfying 
the inequality lul < uj < ujh at all points on its path 
through the plasma can leave the plasma and contribute 
to the wake field spectrum. Within the geometric optics 
approximation, the path of different parts of the wake 
field is determined by the ray equations 
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~dt 



= -d z W , 
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where d/dt = d t + v g d z . Specifically, for a time- 
independent medium, the right hand side of the last equa- 
tion is zero and the wake field propagates, with the local 
wake field group velocity v g = dkW with unchanged fre- 
quency. 

It is natural to first consider the density profile recon- 
struction from a wake field in the case of a quasi-static 
plasma. The proposed technique can be applied not only 
to truly static density profiles but to the large class of 
phenomenons varying on a time scale significantly larger 
than the time of propagation of the wake field through 
the plasma slab. Then, sequential laser pulses gives a 
sequence of wake fields spectra from which a "movie" of 
the density profile can be obtained. 

We assume that the generated wake field spectrum is 
measured immediately outside the plasma boundary. In 
consistence with the geometric optics approximation we 
will treat the measured data as a weakly time dependent 
spectrum with well defined sharp (quasi-monochromatic) 
peaks. The data is not necessarily monochromatic at a 
given time, however. A part of the wake field generated at 
an earlier time may overtake some part of the wake field 
generated at a later time, leading to a spectrum with mul- 
tiple sharp peaks. Generally, we can express the data as a 
set of distinct frequencies measured at different times, in 
which case a sharp curve can be recognized, see Fig. 1(a). 
Due to cut-offs and/or resonances, the curve may be dis- 
continuous, and the algorithm below only works until the 
first jump. 

Before presenting the reconstruction algorithm we need 
to introduce some notations. As mentioned, the mea- 
surement results in a function of time, in general multi- 
valued, describing a path uj(t) as illustrated in Fig. 1(a). 
Let the path be discretized into consecutive points num- 
bered i = 0, 1, ...,n, i = being the first measured point. 
Let also the plasma, having length L, be discretized into 
cells, numbered j — 1,2,...,N , of length Az = L/N 
so that the plasma has the plasma frequency at 
(j — 1) Az < z < jAz. In order to distinguish the two dis- 
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cretizations we will use upper indices for quantities dis- 
cretized according to the spectrum-path discretization, 
e.g. fW, anc i lower indices for quantities discretized ac- 
cording to the plasma discretization, e.g. ■ In par- 
ticular we will let Zrj-j = (j — l)Az whereas z^' is the 
"exact" point in the plasma where the wake field corre- 
sponding to the «:th point in the spectrum was generated. 
Furthermore, it is convenient to introduce Az^ l \ which 
is defined in Fig. 1(b). 

The following algorithm can be used to retrace the 
measured frequency spectrum and thereby reconstruct 
the electron density profile: We assume that the points 

0. 1, 2, i — 1 (grey dots in Fig. 1(b)) have already been 
retraced so that the plasma frequency of the cells j + 

1, j + 2, TV (black) is known and the plasma frequency 
of cell j (grey) has yet to be determined. Consider now 
the retracing of the next spectrum-point i, correspond- 
ing to the frequency wW detected at time t^K Given the 
plasma frequency of the already reconstructed cells, the 
path of spectrum point i is known up to cell j, where 
the plasma frequency and thereby the group velocity is 
yet to be determined. Since spectrum point i at least be- 
longs in the neigborhood of cell j, a first guess - correct 
to lowest order - is that w p (j) = ■ Making this guess, 
the path of spectrum point i can be reconstructed a bit 
further by comparing the detection time, with the 
generation time of this particular wake field, and using 
the ray equations (||). In particular it can be decided 
whether this point belongs to cell j or cell j — 1. If it 
belongs to cell j, the algorithm assigns the calculated 
position to point i and proceeds to spectrum point i + 1. 
If it belongs to j — 1, the cell j is filled with spectrum 
points and a definite plasma frequency is given to this 
cell, before the algorithm continues in the same manner 
with reconstructing cell j — 1 and the path of spectrum 
point i + 1 . For density variations on a faster time-scale 
than c/L, the above algorithm does not work. Prelimi- 
nary work suggests that it can be generalized to a more 
rapidly varying regime, however, but it would require 
multiple laser pulses propagating in the plasma slab si- 
multaneously. 

In order to demonstrate the method we numerically 
calculate a wake field spectrum using the ray equations 
(^) starting from an assumed density profile, see Fig. 2(a). 
In the reconstruction algorithm this spectrum is then 
treated as experimental data. In this case we consider 
a plasma magnetized such that u c = 1.1 x w Pimax , where 
<^p,ms.x is the maximum value of the plasma frequency. 
For simplicity we normalize such that L = 1, c = 1, and 
we let the laser pulse enter the plasma at t = and exit 
at t = 1. The spectrum is obtained by treating the wake 
field as a collection of point particles - jumping into exis- 
tence as the laser pulse propagates through the plasma - 
each moving with the local group velocity and the initial 
wavenumber and frequency determined by Eq. (pp. Since 



the density profile is known, in contrast to the case of re- 
tracing the spectrum, it is straight forward to propagate 
these "particles" by means of the ray equations @. 

Retracing the spectrum according to the algorithm 
presented above results in a density profile that can be 
compared with the one we assumed, see Fig. 2(b). A 
small numerical error - that can be removed with a finer 
discretization - can be seen. Note that the entire plasma 
profile cannot be reconstructed. The left most points in 
Fig. 2(b) are missing. This is because the wake field 
generated in this region of low density cannot propagate 
through the plasma since there is a cut-off prohibiting 
this. The information of this region is already missing in 
the wake field spectrum. 

Existing techniques for unperturbing plasma den- 
sity diagnostics include reflectometry and interferometry 
methods as well as incoherent scattering ||. It is most 
natural to compare our proposal with the two former 
methods: Reflectometry essentially makes use of the re- 
flection of electromagnetic waves against a cut-off occur- 
ring at a point where the frequency of the electromagnetic 
wave equals the local plasma frequency. The plasma fre- 
quency as a function of position can be derived from time- 
of-flight measurements for a suitable spectrum of electro- 
magnetic waves. Applying interferometry, a laser beam 
is split-up into two (or more) beams. One of the beams 
passes through the plasma and experiences a phase shift 
relative to the beam propagating in vacuum. This phase 
shift - obtained by letting the beams interfere - is re- 
lated to the density averaged along the path through the 
plasma. 

In many respects our new method compare favorably 
with existing techniques, provided the basic criterion of 
a strong external magnetic field is fulfilled. Some ad- 
vantages are: 1) Coaxial profile information: Interfer- 
ometry can at best provide information about the den- 
sity profile in a direction perpendicular to the penetrating 
laser beams (if multiple beams are used) whereas the pro- 
posed method reconstructs the density profile along the 
path of the laser pulse. 2) High spatial resolution: The 
requirement that the length scales of inhomogeneities 
must be much larger than the local plasma wave length 
A p = 2itVg/uj p sets the spatial resolution limit to roughly 
10 8 n _1 / 2 m for a strongly magnetized plasma (uj c ~ uj p ), 
where n is the electron number density in units m~ 3 . 3) 
High temporal resolution: The temporal resolution limit 
is equal to the time that separates sequential laser pulses 
probing the dynamics, and thus t rcs ~ AL/c, where AL 
is the spatial separation of the laser pulses. For use of the 
"static" density profile reconstruction scheme AL must 
be comparable to the width of the plasma column, imply- 
ing that i res is not longer than 10 _8 s, for typical labora- 
tory plasma dimensions. 4) Density minimum access In 
our method the density profile can be measured even in 
the vicinity of a local density minimum, in contrast to re- 
flectometric methods for which only monotonous profiles 
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can be fully reconstructed. This feature together with 
points 2) and 3) makes it possible to study the density 
evolution in detail of various dynamical processes, e.g. 
wave propagation and instabilities, using our method. 

Given a suitable plasma - a regime close to thermonu- 
clear fusion conditions seems appropriate - the equip- 
ment necessary to apply the proposed scheme is firstly a 
laser that can produce powerful and short pulses. Com- 
mercial Ti:sapphire lasers typically fulfill these require- 
ments by far. Secondly, the detector need to operate over 
a broad frequency band and produce spectrums that are 
time-resolved on (roughly) the time-scale L^/c, where 
L;h is the density inhomogeneity scale length. Similar 
detector requirements have been fulfilled in the applica- 
tions of ultra-short pulse reflectometry (J], however, and 
thus we deduce that it is not necessary to develop new 
technology for our scheme to work. 

There are a numbers of issues not addressed in our den- 
sity reconstruction scheme. Firstly, what happens if the 
external magnetic field is not homogeneous? This is sim- 
ple provided the magnetic field profile is know, since our 
static reconstruction scheme is essentially unaffected by 
this change, and we only need to include the correspond- 
ing variations in the dispersion relation (^). Secondly, 
what happens if the profile to be reconstructed has two- 
or three-dimensional spatial variations? Unfortunately - 
unless the extra variations occur on longer spatial scales 
- this complicates the application a lot, at least from a 
practical point of view. Two parts of the wake field gen- 
erated at different positions will follow different orbits 
in the plasma and exit at different locations. By using 
multiple lasers and detectors placed at suitable positions, 
it seems possible that 3-D density profiles can be recon- 
structed in principle, but to address this question in any 
detail is beyond the scope of this paper. Finally we point 
out that the question of the usefulness of our scheme is 
unlikely to be settled by purely theoretical arguments, 
and instead we want to encourage experimental investi- 
gations of wake field generation in magnetized plasmas. 
A first interesting step, before complete density profile 
reconstructions are performed, would be to experimen- 
tally verify that the maximum frequency of the wake field 
coincide with the maximum value of the local plasma fre- 
quency. 
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I. FIGURE CAPTIONS 



Fig. 1 Illustration of the discretization scheme, (a) 
Spectrum path discretization (b) Plasma slab discretiza- 
tion 

Fig. 2 Example of a reconstruction, (a) Predicted 
wake field spectrum normalized against the maximum 
frequency (b) . Normalized plasma profiles. The solid 
line marks the assumed profile and the crosses marks the 
reconstructed profile. 
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